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^ ' Abstract 

1^ . Nucleus-nucleus interaction is studied in the framework of the perturbative QCD with 

I Nc ^ oo and a fixed coupling constant. The pomeron tree diagrams are summed by an 

effective field theory. The classical field eqations are solved by iteration procedure, which is 
found convergent in a restricted domain of not too high energies and atomic numbers. The 
', found gluon distributions do not scale, have their maxima close to 2 GeV independent of 

I rapidity and fall towards the central rapidity region. The cross-sections slowly grow with 

' energy due to the contribution from peripheral collisions, where evolution remains linear. 

, Simple variational estimates at higer rapidities confirm this tendency. 

m 
o 

di' 1 Introduction 
(D 

'■^ As discussed in [1] in the perturbative QCD with a large number of colours and a fixed 

coupling constant high-energy nucleus-nucleus interaction is described by the exchange of 
^ ■ an arbitrary number of BFKL pomerons, which interact between themselves via the three- 

pomeron coupling. Corrections due to interactions not reducible to pomeron exchange but 
rather to gluonic exchange are of the order The resulting pomeronic diagrams can 

be classified according to the number of pomeronic loops. Each pomeronic loop gives an 
additional factor So in the high-colour limit only tree diagrams survive. In the case 

of the scattering on the nucleus of a very smal probe (e.g a highly virtual photon) this 
leaves only pomeronic fan diagrams, which can easily be summed to lead to the non-linear 
BFKL evolution equation [2, 3, 4]. This equation, although not soluble analytically, can be 
comparatively easily solved by numerical methods (e.g. [4, 5, 6]). 

For nucleus-nucleus scattering the situation complicates enormously. The basic complica- 
tion comes from the fact that now the pomerons not only split into two but also merge from 
two to one. The tree diagrams now do not reduce to fans but involve other structures, like 
shown in Fig. 1. Still, using methods of the effective non-local field theory, one can sum all 
these diagrams, reducing the problem to the solution of a pair of non-linear field equations 
in the rapidity-transverse-momentum space [1]. Unfortunaltely, contrary to the non-linear 
BFKL equation, these are not evolution equations but rather correspond to a system of 
full-fledged non-linear integral equations, which are very difficult to solve. 

In this paper we make a first attempt at a solution of these equations and try to gain 
some insight into the physical picture of the nucleus-nucleus interaction in this approach. We 
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use two different methods of solution. First we try to find tfie solution by iterative methods 
starting from the fan digrams only. Unfortunately our results show that this method is 
convergent in a very restricted region of not too high rapidities and not too large nuclei. 
Beyond this region the iterations do not converge, indicating some sort of qualitative change 
in the form of the solution (a phase transition?). In relation to this it is worth remembering 
that the primitive Glauber approximation formula for the nucleus-nucleus scattering in the 
tree approximation shows a similar singularity as the nuclei become heavy enough (at A = 64 
for nuclei with a constant profile function within their transverse areas) [7]. To move beyond 
the mentioned limits we tried to use a direct variational method to find the stationary point 
of the effective action, choosing the simplest form of the trial fields. Comparison with the 
exact solution where it can be found by oterations shows that the precision of our variational 
results is not high (of the order ~30%). Still it gives a possibility to see the qualitative 
behaviour of the solutions at very high rapidities. 

Our results show that in the nucleus-nucleus collisions the rise of the effective number of 
gluons becomes still more suppressed than in the non-linear BFKL equation case. In fact 
the variational estimates indicate that it may even go down with the growth of the rapidity. 
However this does not clearly reflect itself on the final nucleus-nucleus cross-sections, which 
continue to slowly rise due to the contribution from the peripheral parts of the nuclei, where, 
due to the the small nuclear density, the cvohition remains practically linear. 

In general the effect of the pomeron interaction on the nuclear cross-section is not very 
impressive. This is a consequence of the fact that the nucleus-nucleus amplitude gets au- 
tomatically unitarized due to cancellations between contributions of different disconnected 
parts. A much greater change can be seen in the contribution of each such part (the eikonal 
function) , which at high rapidity becomes many orders of magnitude smaller than in the pure 
linear BFKL evolution case. 

The paper is organized as follows. In Section 2. we remind our basic formalism to treat the 
nucleus-niiclciis collisions in the perturbative QCD approach, which reduces the problem to 
searching for a stationary point for the action of a certain non-linear and non-local field theory 
In Section 3 we outline our methods to find this stationary and to solve the corresponding 
variational field equations. Section 4 presents our numerical results which are discussed in 
Section 5. 

2 AB-cross-sections and effective field theory 

At fixed overall impact parameter b and (high) rapidity Y the nucleus-A-nucleus-B total 
cross-section is given by 



Here the eikonal function T is a contribution from the connected part and is an integral over 
two impact parameters bA and 6^ of the collision point relative to the centers of the nuclei 
A and B: 



In the perturbative QCD, in the large Nc limit, the eikonal function is given by a sum of all 
connected tree diagrams constructed of BFKL pomerons, which interact between themselves 
via the triple pomeron vertex (with a minus sign). It can be shown that this sum is generated 

by an effective field theory of two fields cf>{y,q) and c/)^{y,q) depending on rapidity y and 
transverse momentum q with an appropriately chosen action S [1] . The action consists of a 




(1) 



T{Y, b) = / d^bAd^bBS^b -bA + bB)T{Y, bA, be) 
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free part Sq, interaction part Sj and external part Se- The free part is given by 

So = 2(</,t|i^(^A + i7^|^) (3) 

where H is the forward BFKL Hamiltonian for the so-cahed semi-amputated amphtudes [8] 
and if is a differential operator in q commuting with H 

K = vyvl (4) 

Symbol (...) means integrating over y and q with weight l/(27r)^ Action 5*0 generates prop- 
agators which are BFKL Green function with operators K^^ attached at their ends. The 
interaction part of the action describes splitting and merging of pomerons: 

5, = ^((^tV0 + ,/.2K,/>t)). (5) 

The coefficient in this term depends on the normalization of the fields. Finally the external 
action is 

SE = -{{wA(t^ + WB<t>^))), (6) 

where wa,b describe the interaction of the pomerons with the projectile and target. If the 
colour distribution in the target is given by 

PA{r)=g''ATA{hA)p{r), (7) 

where p(r) is the colour distribution in the nucleon and Ta is the target nucleus profile 
function, then 

WA{y,q) = 5{y) j (frr^pA{r) = 5{y)wA{q), (8) 

the 5 function indicating that the target is taken to be at zero rapidity. Function wsiy, q) is 
given by a similar formula with S{y) substituted by S{y — Y) wher Y is the rapidity of the 
projectile. 

The classical equations of motion which follow, multiplied by {1/2)K~^ from the left, are 
(^^ + H^4>+ ^(</.2 + 2K-^^^K^) = \k-^wa (9) 

and 

d \ X 2a'iN,. / +2 I +\ 1 1 , ^ 

-— + Hj(t>^ + ^-^[(P^ +2K-^(pK(t>^) = -K~^WB. (10) 

From the J-like dependence on y of the external sources it follows that the equations can be 
taken homogeneous in the interval < y < Y, action of the external sources substituted by 
the boundary conditions 

ct)iO,q) = ^K-'wA{q), cj>\Y,q) = ^K-'wB{q). (11) 

The eikonal function T{Y, bA, bB) is just the action S calculated with the solutions of Eqs. 
(9)and (10), cl)a and <^J,: 

T{Y,bA,bB) = -S{4>cl,4l}- (12) 
Using the equation of motions one can somewhat simplify the expression for S. Indeed 
multiplying the first equation by 2K(p\ the second one by 2K(p, integrating both over y and 
q and summing the results one obtains a relation 

2So + 3Si + Se = 0, (13) 
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which is vahd for the classical action, that is, calculated with the solutions of Eqs. (9) and 
(10). Using this relation we can exclude, say, Sq from (12) to find 

T{Y,bA,bB) = l(Si{<l>ci,<l>li} - SEicPclAll})- (14) 
The dependence on 6^ and comes from the boundary conditions (11). 



3 Methods of solution 

3.1 Final formulas for calculation 

To solve Eqs. (9 and (10) we first rescale the rapidity and fields to pass to variables known 
from studying fan diagrams [3, 4]: 

y^y/a, H ^ aH, '^^2^^' "^^"^202^^' ^-^^^ 

where standardly a = agNc/ir. In these variables, for < y < Y, the equations of motion 
have the same form (9) and (10) without the coefficient before the nonlinear terms and with 
zero right-hand side. All parts of the action aquire a common coefficient 1/(2q!^): 

S/ = 2^((**'-'f* + /ii'*')> (IT) 

and 

Se = -^{cP^KcP^Siy) + d{y - ¥))), (18) 

where we expressed the external sources via the boundary values of (p and (p^ . Note that the 
expression for Sq assumes integration over all values of y, so that the derivative in y generates 
(5-like terms which partially cancel with the external part of the action. If one symmetrizes 
Sq in (j) and (j)^ then these terms cancel exactly one half of 5^. This implies that taking in 
the integration over y in the interval < y <Y one has to take the total action as 

S = So + Sj + ^Se. (19) 

The operator appearing before the second non-linear term in (9) and (10) can be 
represented as an integral operator in the transverse momentum space with a kernel [1] 



(20) 



where 5>(<) = max(min){gi, The operator K contains the 4th derivative in q. To 
simplify it we present it as a product 

K = L^L, L = q^Vl. (21) 

In logarithmic variables L reduces to the 2nd derivative. If 



(22) 
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then 

1 52 

Using this and integrating by part to exclude higher derivatives we can express the compli- 
cated 2nd non-linear term in the equations via finctions (j) and 0^ and their first and second 
derivatives in t. Denoting 

^ = 01, ^ = ^2 (24) 

and similarly for we find 

+ dti(t)2 ((1 - z)4 + 2/3(2z - 1)<^1 - 4p'^z4) } , (25) 

where z = P{t — ti). The conjugated term has the same form with (f) ^ 4 ■ 

Calculating the action one can split K into a pair of operators L acting on factors de- 
pending on (f) and 4 ■ this way one obtains 

(symmetrized in (f> and 4)^ 

Si = < 202(4'/' + 0?) + h.c>, (27) 

Se = < 4h{Siy) + 5iy -Y))>. (28) 

3.2 Boundary conditions 

To fix our boundary conditions we use our experience with the non-linear BFKL equation to 
study the nuclear structure functions [4, 6]. The adequate initial values for ^(y, q) were taken 
there from the Golec-Biernat-Wuesthoff distribution, fitted to the proton data at compara- 
tively low values of x [9], which was duly eikonalized for a nucleus target. In fact cikonaliza- 
tion implies including terms of higher orders in 1/N'c, outisde the precision of the approach. 
Also it is not clear how to generalize eikonalization procedure to the nucleus-nucleus case. 
For both of these reasons our first choice (I) for the initial values is the non-eikonalized 
Golec-Biernat- Wuesthoff distribution for the nucleus: 

«0,,) = -iaEi(-^). (29) 

Here a carries information about the nucleus and impact parameter 

a = aoTAibA), (30) 

(To = 20.8 mb and q is in GeV/c. The value of (/)^(y, g) was taken in the same form with 
Ta(&a) TsihB)- To study a possible influence of the form of the initial distribution in q 
we also used an alternative choice (II) with the same infrared behaviour and point where the 
gluon distribution is peaked but a much slowlier fall of the distribution at large q 

,(0,,) = -l„,„(l + ^). (31) 
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3.3 Iterative solution 

Our first method to find the stationary point of the action has been to solve the classical 
equations of motion iteratively. We have chosen the sum of pure fan diagrams as a starting 
function for iterations. In practice this means that we first solve the equations with the non- 
linear term mixing (j) and (j)^ put to zero. These solutions seve as an input for the iterations 
(f)^""^ and Then we find next iterations from the equations 

^l. + ^{n+l) ^ ^(n+l)2 ^ 2ir- V^^"^i^</'^"^ = (32) 

and 

{^~^ + H^ 0^^"^'^ + <^^^"^'' + 2K-'<I>(^^K<I>^^''^ = (33) 

For each iteration we have only to evolve the initial function from y = Oto y = Y, rather than 
solve the equivalent pair of two dimensional non-linear integral equations, which considerably 
diminishes computer time. 

Unfortunately our calculations show that this method works only for rather small values of 
the participant atomic numbers and rapidity Y. Obviously the maximal value of the factor 
a entering (29) or (30) is achieved at 6 = and for max{^, B}. With the Woods-Saxon 
nuclear density for Pb-Pb, Cu-Cu, Al-Al, 0-0 and C-C collisions this max{a} is found to 
be equal 2.20, 1.53,1.11, 0.88 and 0.83 respectively. Our calculations show that for a < 2.2 
the described iteration procedure is convergent only up to F = 1.1 for choice I of the initial 
distribution and up to y = 1.8. for choice II. For lighter nuclei, taking a < 1.0 we find that 
iterations converge up to y = 1.3 for choice I and up to y = 2 for choice II. The physical 
values of these rapidities depend on the chosen value of a. With a = 0.2 the above numbers 
are to be multiplied by factor 5, but still remain rather low: for, say, 0-0 collisions iterations 
allow to move only up to rapidities of the order 10, that is cm. energies of the order 150 
GeV and for Pb-Pb collision the upper limit for the cm. energy lowers to ~ 90 GeV. 

3.4 Variational solution 

A clear alternative is obviously to try to directly find the stationary point of the action 
choosing some trial fields 4'{y,q) and 4'\y,q) which satisfy the boundary conditions. The 
difficulty of this approach is related to the fact that the action can have more than one 
stationary point. In our first attempt we have chosen the simplest form for the trial fields 
with y and q dependence factorizcd. Moreover for the the y dependence we chose a simple 
exponential one, with a variable slope A, so that our trial fields have the form 

cP{y,q) = e'^ymq), <PHy,q) = e'^^^'-y^^Y^q) (34) 

The boundary functions ^(0, q) and ^^(y, q) were taken according to (29) and (30) for variants 
I and II. The only variational parameter A was chosen to give the minimal value for action S. 
Note, that with the fields having a simple analytic form, the necessary derivative functions 
entering Eqs. (26)-(28) are easily found, so that calulating the action reduces to just doing 
two independent integrations over y and q. With these trial fields the solution for A always 
exists for any values of Y and parameter a in Eqs. (29)- (30) and moreover it corresponds 
to tne minimum of the action. The quality of this approximation can be checked at y and a 
where the exact solution can be found perturbatively. At y = 1 and a = 1 the exact values of 
action S (without factor l/(2a^)) are —0.0120 and —0.0370 for variants I and II respectively. 
The variational values obtained with (34) for these two variants are —0.0100 and —0.0262. 
As one observes the precision is not very high (especially for variant II). Still we hope that 
the variational approach might give some indication about the behaviour of the cross-section 
and eikonal functions at large values of Y at which we cannot obtain the exact solution. 
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4 Numerical results 

We first report on the iterational solution of the field equations (32), (33), which as mentioned 
is convergent at not too high values of y, ^ and B. We chose to study 0-0 scattering 
{A = B = 16) using choice II of the initial functions, which allowed us to obtain the solution 
up to y = 2. Presenting our results we first consider the gluonic density, which can be related 
to functions Lcj){y,q) = h{y,q) and L(f)\y,q) = h'^{y,q). Indeed as follows from the study of 
the non-linear BFKL equation the gluon density of a single heavy nucleus is given by [4] 

dxG{x,q) _ 1 

We conjecture that the gluon density in the nucleus-nucleus collision at rapidity y will be 
given by a similar formula with contributions from both nuclei. For central collisions then 



In Fig. 2 we present our solution for h{y,q) at y = 2, 6^ = 6^ = at different stages 
of evolution: y = 0, 0.5, 1.0, 1.5 and 2. For comparison we show in Fig. 3 the same 
function which is found from the non-linear evolution equation for a single nucleus (only fan 
diagrams). The difference between Fig. 2 and 3 comes from the influence of another nucleus 
on the evolution process. As one can conclude, this influence is quite strong. Whereas the fan- 
diagram density steadily shifts towards higher momenta more or less preserving in its shape, 
the density for the nucleus-nucleus collision practically does not move until y = 1, its peak 
dramatically falling and its low momentum tail visibly growing. Only as late as at y = 1.5 
one notices some slow shift towards higher momenta, which becomes more pronounced at 
y = 2. Still at y = 2 its peak lies at q = Qs = 2 GeV/c ("saturation momentum"), whereas 
for single nucleus it is found at = 8 GeV/c. In Fig. 4 we illustrate the total gluon density 
in 0-0 collisions at y = 2 given by the sum h{y, q) + h{Y — q) up to a factor depending on the 
coupling constant value. This density has its peak at the point close to 1 GeV/c practically 
independent of rapidity. The height of the peak is falling towards the central region. Some 
diffusion towards small momenta is observed. It can however hardly be compared to the 
diffusion for the pure BFKL evolution illustrated in Fig. 5 for the same initial function and 
same region of y. 

A clear physical observable is of course the total nucleus-nucleus cross-section obtained 
by the integration of (1) over all impact parameters. We show it for 0-0 scattering at y < 2 
in Fig. 6. To compare we present also the cross-sections corresponding to a single BFKL 
exchange. These latter are naturally larger but the difference is not at all dramatic, reaching 
some 18% at the maximal rapidity y = 2. This is understandable, having in mind that Eq. 
(1) actualy automatically unitarizes the amplitude and leads to very similar results even for 
very different eikonal functions provided they are large. This can be clearly seen from the 
comparison of the eikonal functions at 6 = in Fig. 7. At y = 2 the pomeron interaction 
in the nucleus-nucleus collisions reduces it by a order of magnitude, although it still remains 
large, ^ 200. Some spreading of the gluon distribution into the low momenta domain visible 
in Fig. 4 in the central rapidity region makes one think that the results may be rather 
sensitive to the infrared region and so strongly dependent on the infrared cutoff. Such a 
dependence indeed exists but is not so strong. With the infrared cutoff at kmin = 0-3 GeV 
we obtain at y = 2 ao-o = 4.18 bn and To_o(0) = 179 whereas without cutoff we have 
(To-o = 5.36 bn and To_o(0) = 208. As one observes the cross-section results more sensitive 




(36) 



Note that in the considered symmetric case </>^(y, q) = (f){Y — y, q 



and so h^{y, q) = h{Y — y, q). 
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to the infrared cutoff, which is a result of peripheral collisions, where evolution follows the 
linear BFKL equation. 

Our variational results for both variants of the choice of the initial functions are presented 
in Figs. 8-11. Since in this approximation we are not restricted to small values of A and 
B, we show our results for Pb-Pb collisions {A = B = 207) In Figs. 8 and 9 we show the 
total cross-sections for variants I and II of the initial fields. They are again compared to the 
cross-sections corresponding to the single BFKL exchange. All the cross-sections steadily rise 
with Y. However this rise seem to be very weak for the choice I of the initial function. The 
single BFKL exchange naturally leads to larger cross-sections and the ratio of these to the 
cross-sections with pomeronic interaction rises with Y reaching values 2 and 1.5 at F = 6 for 
variants I and II respectively. However this difference is far larger for the eikonal functions, 
shown in Figs. 10 and 11 at 6 = 0. The eikonal function for a single BFKL exchange rises 
up to values of the order 10^ at Y = 6, whereas with pomeronic interactions we find values 
around 100 or 1000 for variants I and II. It is remarkable that, with the pomeronic interaction 
switched on, the eikonal function actually diminishes with y for central collisions. Therefore 
the rise of the cross-section is totally due to peripheral collisions, where, with a low nuclear 
density, the non-linear effects are small and the fields grow according to the pure linear BFKL 
equation. 

5 Discusssion 

We have made the first attempt to solve the equations which describe nucleus-nucleus scat- 
tering in the framework of the perturbative QCD with a large nucmber of colours and a fixed 
coupling constant. The natural iterative approach has been found to converge in a restricted 
domain of not too high scaled rapidities and atomic number of participants. Physical ra- 
pidities covered by the convergence range depend on the value of the coupling constant. For 
as = 0.2 they are not greater than 10 for 0-0 collisons and not greater than 9 for Pb-Pb col- 
lisions. The solutions in this range of Y has been found to generate the gluon density which 
falls with y (towards the central rapidity region) and somewhat spreads into the infrared 
region of transverse momenta, its maximum staying around 1 GeV/c. It radically differs 
from the density of the isolated nucleus, which is known to steadily shift towards higher 
momenta, the hight of its peak practically independent of y. Both the eikonal function and 
the total cross-section are found to be damped as compared to the single pomeron exchange 
(by an order of magnitude for the eikonal function). The latter has been found to actually 
fall with energy for central collisions. However the total cross-sections rise with energy due 
to peripheral collisions where non-linear effects are naturally small. 

Unfortunately we have not been able to find the solutions outside the mentioned restricted 
domain of rapidities and atomic numbers. We do not know what sort of singularity occurs 
at the boundaries of this domain and even if the solutions of our equations exist at all. It is 
possible that a sort of phase transition occurs at this boundaries, so that the equations have 
to be changed. 

Just to see some qualitative features of a possible solution at high rapidities and atomic 
numbers we applied a simple variational procedure, approximating the fields by certain simple 
trial functions with a single parameter to be determined from the stationary point equation. 
It is hopeless to expect to study the gluon distribution from such a simple approach. One 
expects more reasonable answers for the eikonal and especially for the total cross-sections 
which are weakly dependent on moderate variations of the fields. Our results seem to indicate 
that with the further rise of energy the cross-sections continue to grow slowly, much slowlier 
than with a single BFKL exchange. The eikonal function in the center continues to fall, very 
slowly in variant II for the initial function and rather fast for variant I. 
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The main lesson to be learned from these first calculations is that the dynamics of nucleus- 
nucleus collisions is much more complicated than for collisions of a small probe on a single 
nucleus. The gluon densities we have found have a much more complicated form than in 
the latter case when they scale with the saturation momentum which grows with energy as a 
power. No scaling of this sort have been observed. The remaing problem is to understand the 
reason of the breakdown of the iterative solution at a certain value of energy and/or atomic 
number and to try to move beyond this value. 
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7 Figure captions 

Fig. 1. Some pomeron tree diagrams summed by our effective field theory. 

Fig. 2. The gluon distribution corresponding to the field (p{y, q) with initial function II. 
Numbers show rapidities y = 0, 0.5, 1., 1.5 and 2. 

Fig. 3. Same as Fig. 2 for the field evolving according to the non-linear BFKL equation 
(only fan diagrams). Curves from left to right correspond to y = 0, 0.5, 1. 1.5 and 2. 

Fig. 4. The total gluon density in 0-0 collisions. Numbers show rapidity distances from 
the target or projectile: y = 0, 0.5 and 1. 

Fig. 5. Same as Fig. 2 for the field evolving according to the linear BFKL equation. 
Curves from bottom to top correspond to rapidities y = 0, 0.5, 1., 1.5 and 2. 

Fig. 6. The total cross-section for 0-0 collisions with initial conditions II (lower curve). 
The upper curve corresponds to a single BFKL exchange. 

Fig. 7. The eikonal function at 6 = for 0-0 collisions with initial conditions II (lower 
curve). The upper curve corresponds to a single BFKL exchange. 
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Fig. 8. Variational estimates for the total cross-sections for Pb-Pb collisions at high 
rapidities with initial conditions I (lower curve). The upper curve corresponds to a single 
BFKL exchange. 

Fig. 9. Same as in Fig. 9 for initial conditions II. 

Fig. 10. Variational estimates for the eikonal function at 6 = for Pb-Pb collisions at 
high rapidities with initial conditions I (lower curve) . The upper curve corresponds to a single 
BFKL exchange. 

Fig. 11. Same as in Fig. 10 for initial conditions II. 
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